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ABSTRACT 

A complete model of helium-like line and continuum emission has been incor- 
porated into the plasma simulation code Cloudy. All elements between He and 
Zn are treated, any number of levels can be considered, and radiative and coUi- 
sional processes are included. This includes photoionization from all levels, line 
transfer including continuum pumping and destruction by background opacities, 
scattering, and collisional processes. The model is calculated self-consistently 
along with the ionization and thermal structure of the surrounding nebula. The 
result is a complete line and continuum spectrum of the plasma. Here we focus 
on the ions of the He I sequence and reconsider the standard helium-like X-ray 
diagnostics. We first consider semi-analytical predictions and compare these with 
previous work in the low-density, optically-thin limit. We then perform numerical 
calculations of helium-like X-ray emission (such as is observed in some regions of 
Seyferts) and predict line ratios as a function of ionizing flux, hydrogen density, 
and column density. In particular, we demonstrate that, in photoionized plasmas, 
the i?-ratio, a density indicator in a collisional plasma, depends on the ionization 
fraction and is strongly affected by optical depth for large column densities. We 
also introduce the notion that the i?-ratio is a measure of the incident continuum 
at UV wavelengths. The G-ratio, which is temperature-sensitive in a collisional 
plasma, is also discussed, and shown to be strongly affected by continuum pump- 
ing and optical depth as well. These distinguish a photoionized plasma from the 
more commonly studied collisional case. 

Subject headings: X-rays: galaxies — methods: numerical — atomic processes — 
plasmas 



1. Introduction 



X-ray emission lines from high-excitation ions offer a different view of Active Galactic 
Nuclei (AGN) emission line regions than do the strong UV and optical lines that are usually 
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studied. While emission lines such as H I Lya, He I A5876, and C IV A1549 come from gas 
with electron temperature Tg ~ lO^K, the emission lines which are detected in X-ray spectra 
can come from gas that is much hotter and more highly ionized, and which presumably lies 
closer in to the ionizing continuum source. Such X-ray emission lines are often discussed 
in connection with observations of X-ray and UV warm absorbers (eg. Kraemer et al. 
2005; Netzer et al.2003; George et al. 1995). Emission lines from helium-like ions are 
important in this regard, because the ratios of their X-ray lines provide temperature and 
density diagnostics that can be measured from intermediate-resolution X-ray spectra. Several 
of these lines have now been detected in a number of nearby AGN, including Mrk 279 
(Costantini et al. 2005; Kaastra et al. 2004), NGC 5548 (Steenbrugge et al. 2005), NGC 
3783 (Netzer et al. 2003; Kaspi et al. 2002), NGC 4051 (Collinge et al. 2001), NGC 7469 
(Blustin et al. 2003), and NGC 4151 (Kraemer et al. 2005). 

The X-ray forbidden (f) and intercombination (i) lines connect the triplet states to 
the singlet ground state, while the resonance (r) line is due to transitions into the ground 
state from 2p^P. (See the schematic energy-level diagram in Figure [1]). Diagnostic line 
ratios for obtaining temperature, G(Tc) = (/f + /i)//r, and density, R{ne) = U/h, from 
particular X-ray line ratios in plasmas have been discussed by a number of authors (Bautista 
& Kallman 2000; Porquet & Dubau 2000; Porquet et al. 2001; Pradhan & ShuU 1981; 
Pradhan 1985; Blumenthal, Drake, & Tucker 1972; and many others). These were introduced 
by Gabriel & Jordan (1969, 1973) and were originally developed for coronal plasmas, not 
where photoionization-recombination is dominant. Recently, authors have discussed the 
influence of radiative transfer effects on the standard X-ray diagnostics (see Coupe et al. 
2004, Godet et al. 2004, and Bianchi & Matt 2002). 

The theoretical models of the He I isoelectronic sequence presented here have been 
developed as part of the plasma simulation code Cloudy (last described by Ferland et al. 
1998). The model helium atom has been discussed by Porter et al. (2005, 2007) and 
Bauman et al. (2005). Where the sources and methods used in the model ions differ from 
those previously discussed sources and methods of the model helium atom, those differences 
are given in Appendix \^ 

Here we consider the standard X-ray diagnostics from both semi-analytical and numer- 
ical perspectives. In section [2], we derive semi-analytical expressions of line ratios in the 
limit where the lines are formed by recombination and calculate the line ratios for helium- 
like oxygen and iron at appropriate temperatures. We compare our results with those of 
other workers and discuss how radiative transfer and collisions will affect the emission. In 
section [3], we include all relevant physics in a numerical model to calculate theoretical line 
ratios as functions of column density. This illustrates the effect of increasing optical depth. 
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We also calculate helium-like emission as a function of ionizing flux and hydrogen density. 
We present the predicted line ratios for oxygen and discuss how our results affect the in- 
terpretaion of existing and future measurements of the emission lines from AGN. We also 
discuss the effects of continuum pumping. We conclude in section HI 

2. Semi-analytical Calculations 

Gabriel & Jordan (1969, 1973) showed how ratios of intensities of various lines involving 
transitions between the n = 2 and n = 1 levels of helium-like ions could be used to measure 
temperature and density in a coUisionally-ionized plasma. In such a gas the electron kinetic 
energy is roughly equal to the ionization potentials of the species that are present. Collisional 
excitation of the n = 2 levels from ground is the main process that produces the X-ray lines. 

These results do not carry over to the photoionized plasmas considered here, however. 
In this second case the ionization and kinetic temperature are controlled by the radiation 
field. The kinetic energy is much less than the ionization potentials of the dominant stages of 
ionization and the n = 2 levels are mainly populated by recombination from the next higher 
ionization stage. As a result the R and G ratios are not density and temperature indicators, 
but, as we show here, are determined by the ionization and column density. 

In the following we will consider helium-like oxygen (recombination onto O^^ forming 
O^^) and iron (recombination onto Fe^^+ forming Fe^'^+) in detail. Tests show that, for an 
AGN continuum and solar abundances, the kinetic temperature is roughly 500,000 K and 
10'' K for these two ions. We will assume these temperatures in the estimates that follow. 

2.1. Level Populations 

The theoretical intensities of our lines of interest (transitions from n = 2 to the ground 
state) depend on the populations of the upper levels of each transition: 2p ^P, 2s ^S, and 2p ^P. 
In a low-density, pure recombination plasma, the time-steady balance equations for the pop- 
ulations of individual levels involve only effective recombination coefficients and transition 
probabilities. We consider that limit here. We define a'/^^Z^T) as the effective recombina- 
tion coefficient (including dielectronic recombination) into level i of the He-like ion of nuclear 
charge Z at temperature T. We calculate a1^^ {Z, T) with the use of cascade probabilities, as 
outlined by Robbins (1968). (See also Osterbrock & Ferland 2006, page 85 for a discussion 
of effective recombination.) We include nLS*- resolved terms with principal quantum number 
n < 20 and a series of 'collapsed' n-resolved levels with 20 < n < 100 (as discussed in Porter 
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et al. 2005 and Bauman et al. 2005). See Appendix|X]for a discussion of our atomic sources, 
and Figured] for an energy level diagram. Note that in Figured] and in the text, we refer to 
the UV lines 2p^Pj — 2s and 2p — 2s as UV3j and UVl, respectively. The individual 
components of the intercombination doublet 2p^Pj — Is are designated by ij. 

Ions can radiatively decay from the 2p ^P term via the electric-dipole transition 2p ^P — 
2s and via the intercombination doublet 2p^Pj — Is^S. The intercombination transition 
probability, A2pSp-^^is^s{Z)y (approximately oc Z^) increases with increasing nuclear charge, 
Z, faster than the electric-dipole transition probability, A2p3p_2s3s{Z), (approximately oc Z"^) 
so that, for Z > 6, the intercombination line is stronger than the electric-dipole transition. 
The transition probabilities to ground from the separate 2p^Pj levels are vastly different. 
The transition 2p^Pq — Is is strictly forbidden (for one-photon transitions by the selection 
rule forbidding J = -v^ 0, and for two-photon transitions by the rule requiring AL = 0, 2) 
while the transition probability A2p3p-^^-is^siZ) is the fastest, about a thousand times faster 
than A2p3p2-isis{Z) for carbon and nearly 7000 times faster for iron. These large differences 
in transition probabilities make it necessary to solve for the populations of the separate 
j-levels separately. The transition probabilities A2p3p._2s3s{Z) also depend upon j. These 
differences are only as large as about 25% for carbon, while in the case of iron the transition 
probability for j = 2 is 3-4 times the j = and 1 values. The populations, n2p3p^, of the 
2p ^Pj term^ in the low density limit are as follows: 

n2p3p,{Z,T) _ go allip^{Z,T) ^^^^^ 



ne nz+ A2p3p^_2s3s{Z) 
n2p3p,iZ,T) _ g,all{p^{Z,T) 



^11+ [A2p3p^-2s^s{.Z) + A2p3p^_isis{.Z)\ 

n2p-3p,iZ,T) ^ 92 <{p„(^,r) 

ne nz+ [A2p3P2-2sSs{Z) + A2p3P2-ls^s{Z)] 



[cm^] (2) 
[cm% (3) 



where nz+ is the density of the hydrogen-like ion, T is the electron temperature, is the 
electron density, and we have used a'^p3pXZ,T) = gj al^ip^{Z,T) and A2pfPo-isis(^) = 0. 

Ions can radiatively decay from the 2p ^P level to either the ground state, 2s ^S, or 
2s '^5'. The decay to ground dominates for all Z but is highly susceptible to optical depth 
effects since it is a resonance line - this is discussed more below. The rate of the transition 
to 2s '^5' (approximately oc Z'') increases with increasing Z faster than the decay to 2s 



^The S'T-mixing discussed in Bauman et al. 2005 for helium becomes important for the high Z ions of 
the sequence as the physical 2p^Pi level is mixed roughly 30% with the LS 2p^P level in the case of iron. 
We do not consider this effect here. See Lin et al. (1977). 
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(approximately oc Z"^) so that the two transition probabihties are comparable in the case of 
helium-like iron. The population, n2pip, of the 2p^P term is 

He nz+ [A2p^P~ls^s{Z) + A.2plp-2s^s{Z) + ^2p iP-2s ^5 (^)] 

Ions can radiatively decay from the metastable 2s level to the ground state via two- 
photon and magnetic-dipole transitions. The magnetic-dipole decay dominates for all Z, 
and we accordingly ignore the two-photon decay in the semi-analytical discussion. The 
population, 712^35, of the 25^5" level is given by 

Equations [T]l5] are valid when triplet-singlet exchange collisions (which change the spin 
of one of the electrons) can be neglected. To first order, exchange collisions add an additional 
term in each equation. Equation [5l for example, would be modified as follows: 

n2s3s{Z,T) al%{Z,T) 3 

= / ' A 7^7^ [cm^]- (6) 

He nz+ A2s3s^ls^S{Z) + 2^ q2s-iS-nl^L[Z,T) 

n,L 

where q2s ^s-ni^L{Z,T) is the collision rate coefficient (in units cm'^ s^^) from 25^5* to the 
singlet term nl^L. We define qtot{Z,T) = J2 123^3 -ni^L{Z,T). Colhsions are negligible if 

n,L 

Tie qtot{Z,T) is small relative to A2s3s-is'^s{Z). Considering oxygen at 500,000 K, a temper- 
ature near where the hydrogen-like ionization state peaks, we find ^2s35_isi5(6) = 10^ s~^ 
and gtoi(8, 10^'^ K) 4 x 10~^° cm^ s~^, corresponding to a critical electron density of 
2 X 10^^ cm^. For iron at a temperature of 10'' K, we find a critical density of 2 x 10^^ cm'^. 

Because the transition probabilities (approximately oc Z^^) increase with increasing Z, 
while the collisional rate coefficients decrease with increasing transition energies (approxi- 
mately oc Z^), the critical density will be even larger for heavier ions of the sequence. We 
restrict our semi-analytical calculations to lesser densities and neglect exchange collisions. 



2.2. Emissivities 

The total emission. An jx/ue nz+, of a recombination line with wavelength A is 

Ue nz+ A He nz+ 
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where Aui{Z) is the transition probabihty of the transition and Z"*" is the parent (hydrogenic) 
ion. The dimensionless ratios of the forbidden and resonance emissivities to the intercombi- 
nation emissivity follows: 



R 



JiiZ,T) 
Ji{Z,T) 



Ai 
Af 
Ai 
A. 



A2s^s^is^s{Z) n2s3s{Z,T) ^ A2s3s-is^s{Z) n2sSs{Z, T) 
_A.2p-iPi-is^siZ) n2p3p^{Z,T) A2p'ip2-is^siZ) n2p-Ap^{Z ,T) 
^2pip-isi5(^) n2pip{Z,T) ^ A2pip-isis{Z) n2pip{Z,T) 



(9) 



_A2p'-ip^-is^s{Z) n2p'ip^{Z,T) A2p3p2-isT^siZ) n2p3.p^{Z,T) 
where Ai, Af, and Ar are the wavelengths of the intercombination, forbidden, and resonance 
transitions, respectively and we have given the label "L" to the ratio of resonance to inter- 
combination emission. The G-ratio, defined as the sum of the intercombination and forbidden 
emission divided by the resonance emission, is given by 

jf(Z,T)+ji(Z,T) R+l 



G 



UZ,T) 



(10) 



2.3. Comparison with Previous Work 

The results we obtain from the above analyses are in good agreement with previous 
work. At a temperature of 500,000 K, typical of a plasma in which recombination onto 
helium-like oxygen is important, we find from Equation [8] the O VII ratio R = ji/ji = 4.2. 
At the same temperature. Figure 8 of Porquet & Dubau (2000) indicates R ~ 4.3 in the 
low-density hmit, while Figure 3 of Bautista & Kallman (2000) suggests R ~ 4.2. The value 
calculated here agrees very well with those two values. We find from Equation (TD] the ratio 
G = (jf + ji)/jr = 5.0. Porquet & Dubau find (as taken from their Figure 7) G ~ 4.8 when 
the hydrogen-like ionization stage of oxygen is at its peak. Figure 4 of Bautista & Kallman 
reports G ~ 5.0 at the conditions considered here. Again, the present value agrees very well 
with values found by other workers. 

For the case of iron, at a temperature of 10'' K, we find G = 4.2 and R = 0.68. The 
results from the Bautista & Kallman (2000) work are G = 4.7 and R = 0.71 (Bautista, 
private communication). These results are also in good agreement with our results. 



2.4. Optical Depth Effects 

Introducing escape probabilities, which modify transition probabilities to yield effective 
transition probabilities, the emission per unit volume is given by 

he 

47r /a = y n^iZ,T) A^Z) [ergs cm'^ s-^]. (11) 
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(See Elitzur 1992, for a discussion of escape probabilities.) Small optical depth corresponds 
to an e of unity, and e ~ r^^ when r is large. The optical depth of a line increases as the 
column density of the ion increases but also depends upon the wavelength and transition 
probability of the line. The effect of optical depth on line ratios is dependent upon which 
line becomes optically thick first (or which escape probability becomes less than unity first) 
as column density increases. It is important to note, however, that, in addition to the 
explicit escape probability dependence in Equation [TTl there is also an implicit dependence 
in the population n^j(Z, T)|§. For an upper level with only one significant decay mode, these 
escape probabilities cancel in the low density limit (see, for example, n2s3s{Z,T) given in 
Equation [5]) , and the line is said to be effectively optically thin. For an upper level with more 
than one significant decay mode, the escape probabilities do not cancel out, so line emission 
from these levels is affected by optical depth. For example, the ratio of the forbidden to ii 
intensities (equal to the i?-ratio where the 12 line can be neglected, as with O VII) varies as 
follows 

f_ A2p3p^_isisiZ) ej^ + A2p3p^_2s3s{Z) euv3^ 

where and euvsi are the escape probabilities of the ii and UV3i lines, respectively. See 
Appendix [B] for a demonstration that the intercombination line can become optically thick. 

The intercombination rate dominates in Equation [12] for all Z > 6. The intercombination 
line also becomes optically thick well before the UV3i line. (Optical depth in the UV3j lines 
is also discussed in Appendix [B]) For oxygen at 500,000 K, the i?-ratio begins to be affected 
by optical depth at about A^(0^"^) = 10^^ cm~^. When this happens, as we will see below, 
R will surpass the classical, canonical value as the denominator of the ratio is surpressed. 
For iron, the intercombination lines will not be affected by optical depth at all (this is also 
discussed below). The fact that the intercombination line can become optically thick was 
mentioned by Godet et al. (2004) but was not demonstrated or indicated in any results. 

Recent observational results on the G and i?-ratios in NGC 4151 can be found in Ar- 
mentrout et al. (2007). Kinkhabwala et al. (2002) presented large observed values of the 
i?-ratio (in NGC 1068). The present work demonstrates that optical thickness in the inter- 
combination line is a possible explanation. 

While the G-ratio can depend strongly on optical depth, it is also strongly affected by 
direct continuum pumping of the resonance line. The transition with increasing column 
density from the pumped case to the optically thick case is not well treated with an escape 



^Actually, effective recombinations can also be affected by line optical depths. This is included in the 
numerical calculations below but is a complication beyond the scope of these semi-analytical calculations. 
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probability formalism. Accordingly, these effects on the G-ratio will be neglected in the 
semi-analytical analysis, although they are included in the numerical solutions. 



2.5. Individual Components of the Intercombination Doublet 



Here we present a method by which the optical depth of intercombination lines, and so 
the column density of its ion, can be derived from the relative intensities of the individual 
components of the intercombination doublet. 

The individual components ij of the intercombination doublet are related by the follow- 
ing expression: 



^2p3pi-2g3s(Z)euv3^ -I- y42p3p^_isig(Z)ei^ 

(Z) 



(13) 



where euvsj and euvi are the escape probabilities for the UV transitions 2p ^Pj — 2s and 
2p — 2s ^S, respectively, and ej^ is the escape probability for the ii line. We have not 
included escape probabilities for the 12 line because that line will always have an optical 
depth roughly three orders of magnitude less than the optical depth of the ii line. The 
rightmost factor in Equation [13] cannot be simplified for a general Z, because in both the 
numerator and the denominator the dominant transition probability is a function of Z. 
Again, for the case of oxygen, we have 



4(8, T) _ 21.81 5 3.5 x 10^ 
4 (8, T) ~ 21.8 3 5.5 x lO^e^ 



^.1 X 10^euv3, +5.5 X 10% 
8.4 X lO^euvsa + 3.5 x 10^ 



(0.0011/ei,: 



j.l X 10%v3^ + 5.5 X 10% 
8.4 X lO^euvaa + 3.5 x 10^ 



(14) 



If optical depth effects are not important, li^/Ii^ ~ 0.0080. If optical depths are significant 
(but still small enough that the 3.5 x 10^ in the denominator in the rightmost factor can be 
neglected), we have 



4(8,T) _ O.OOlleuvsi 0.0069 



(15) 



and the ratio will increase with increasing optical depth. If optical depth in the UV3j lines 
is not important (see Appendix [B]), the escape probability (inversely proportional to the 
column density) in the intercombination line can be derived from the individual components 
of the intercombination doublet as follows: 



eil 



0.0011 



4/4 - 0.0069 



f{N). 



(16) 
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A similar but different expression holds for other ions of the sequence. For iron, the ii rate 
is so dominant that it cancels out completely in Equation [131 and the ratio is 

4(26, T) ^ 5 
4(26, T) ^ 3 

The ratio therefore ranges only slightly, varying from about 4/3 for optically thin media to 
5/3 for optically thick media. With sufficiently high signal-to- noise observations, the optical 
depth of the UV32 line could be derived from X-ray observations of the intercombination 
lines. 

It is important to note that this analysis is largely independent of density because 
collisions will only begin to affect the relative populations of the separate 2p^Pj levels for 
densities greater than Ue = 10^^ cm~^ for oxygen (at 500, 000 K) and Ue = 10^^ cm~'^ for 
iron (at 10^ K). 

Resolving the individual components of the intercombination doublet may be difficult. 
The line separation of the doublet corresponds to a velocity range from only about 20 km s~^ 
for carbon to almost 700 km s~^ for iron. If the intrinsic line broadening is small, then the 
intercombination line can be detected in systems with significantly less broadening. The 
individual components could be resolved with a spectrometer having resolving power of at 
least 12000 for carbon and 400 for iron. Current and planned X-ray satellite observatories 
might be able to resolve the individual components for ions as light as magnesium and almost 
certainly for ions as heavy as iron. Unfortunately, optical depth in the intercombination line 
is less likely for heavy ions because they are generally less abundant, and in practice the 
method outlined here may be untenable. The best candidates for resolving the individual 
components of an optically thick intercombination doublet are probably Mg XI, Si XIII, and 
S XV. 



1 

0.24 euv32 + 1 



(17) 



3. Numerical Calculations 
3.1. Summary 

Having established our semi- analytical results in the low-density, optically-thin limit, we 
now consider numerical calculations. In this section, we present the results of Cloudy models 
of helium-like X-ray emission for a wide range of physical conditions. These models span 
the range of emitters that may be found in high-ionization regions of AGN. The emergent 
spectrum is calculated self-consistently with the ionization and thermal structure of the 
line-emitting region. This includes coUisional processes, line transfer, background opacities, 
continuum pumping, and dielectronic recombination. 
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We first consider line ratios for a particular model as a function of hydrogen column 
density. Then we vary the ionizing flux and volume density of the original model and present 
line ratios as contour plots. Finally, we illustrate the effects of continuum pumping by adding 
an additional incident continuum source. 



3.2. Line Ratios as a Function of Column Density 

We calculate a constant density plane-parallel slab with the Korista et al. (1997) ionizing 
continuum incident with hydrogen-ionizing flux $(H) = 10^*^ photons cm~^ s~^ and hydrogen 
density tt-h = 10^ cm~^, corresponding to log f/(H) ^ 0.5. These parameters were chosen 
because helium-like oxygen emission peaks near these conditions for a small column density. 
Solar abundances are assumed (Grevesse & Sauval 1998; Allende Prieto et al. 2001, 2002; 
Holweger 2001). 

In Figure m we plot the O VII R, L, and G ratios as a function of column density. The 
i?-ratio increases by a factor of 2 at large column densities. This is due to an increase in 
the optical depth of the intercombination line (see Equation [12]) and suggests that caution 
should be used in deriving electron density from the i?-ratio in conditions where the column 
density is large and the intercombination line may be optically thick. For low column den- 
sities, R agrees well with the semi-analytical value calculated above. The L-ratio (/i.//i), 
however, falls by a factor of nearly 70 over the same range of column densities. The reason 
is that the resonance line becomes optically thick much faster than the intercombination 
line. The ratio /r//i is ~ 23 times the semi-analytical value at low column densities. This 
is because continuum pumping is dramatically enhancing the resonance line relative to the 
semi-analytical value. For small column densities, photoexcitation is important and the 
plasma is in Case-C conditions (see Ferland 1999 and Baker et al. 1938). For large column 
densities, the plasma exhibits Case-B behavior (Baker & Menzel 1938). The temperature- 
sensitive ratio G = (If + /i)//r varies by more than two orders of magnitude when the column 
density increases from A^h = 10^'' cm~^ to A^h = 10^^ cm~^. This trend suggests that the 
column density must be constrained before using G as an ionization indicator. In the low 
column density limit, the G-ratio is smaller than the semi- analytical value by a factor of 
^ 25, a finding consistent with the resonance line being enhanced by continuum pumping, 
as discussed above. 
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3.3. Line Ratios as a Function of Ionizing Flux and Volume Density 

Next, we calculate a grid of simulations using the same continuum shape and com- 
position but varying both the flux of hydrogen-ionizing photons and the hydrogen volume 
density. We consider two hydrogen column densities, iVn = 10^^ cm~^ and 10^^ cm~^ to show 
this dependence. In Figures HE we plot a number of O Vll intensity ratios as a function 
of hydrogen density and ionizing flux. Note that in a triangular region in the bottom right 
corner of each contour plot, to the right of lognn = 6 and below log<l'(H) = 24, oxygen is 
not ionized enough to produce signiflcant helium-like emission. 

In Figure [3], we plot the ionization fractions of O^''' and 0''+. The plots indicate that in 
the upper left corner of the plotted parameter space both stages are neglible (because oxygen 
is almost entirely stripped of electrons). Moving down and toward the right, the hydrogenic 
ionization state begins to dominate, followed by the helium-like state. In the bottom right 
corner, as noted above, oxygen is too neutral to produce signiflcant helium-like emission. 

In Figure m we plot the ratio R = U/h for both column densities mentioned above. 
Note that, in both panels, for a given ratio and ionizing flux, there is not a unique density. 
Values of R span several orders of magnitude for the entire density range plotted. This flgure, 
combined with Figure [2], demonstrates that the i?-ratio is not a simple density diagnostic. 
Two constant ionization parameter lines are overplotted in each panel in Figure |H One can 
clearly see that, for a given value of i?, a change in ionization parameter leads to an almost 
identical change in the derived density. Considering the i?-ratio as a function of ionization 
parameter causes the appearance of a density dependence in the ratio, but that dependence 
is almost entirely due to the deflnition of the ionization parameter and not to the physics 
involved in the i?-ratio. The i?-ratio is roughly constant at low flux, inversely proportional 
to the flux for several orders of magnitude of the flux, and then constant at very high fluxes. 
In the large column density case (right panel), the ratio is roughly a factor of two greater 
than the corresponding values in the small column density case in the lower left corner of 
the parameter space. Values of the O VII i?-ratio larger than 4 have been observed but were 
previously unexplained by theory (see section [231) • 

The ratio G = {!{ + /i)//r, a temperature indicator in a collisional plasma, is plotted in 
Figure [51 In a photoionized plasma, the ratio is a function of ionization fraction (shown in 
Figure [31) and column density. The ratio depends on both ionization parameter and column 
density. This fact can be seen in the right panel of Figure [3 where the column density is 
large and the G-ratio differs greatly from the small column density case in the left panel. 
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3.4. Effects of Continuum Pumping 

The results presented above are strongly dependent upon continuum pumping. This 
is why the gas reaches Case-C conditions at small column densities (see Figure [2]). Here 
we further illustrate the effects of continuum pumping by considering a second continuum 
component. 

Some of the above results at high flux remain counterintuitive, in particular that the 
i?-ratio in Figure H] should trace ionizing flux (i.e., is parallel to the density axis). This is 
caused by direct continuum pumping of the UV3j lines. A simple calculation illustrates the 
effect. 

We add a 10,000 K blackbody to the original Korista et al. (1997) continuum used 
above (again with hydrogen-ionizing flux $(H) = 10^^ photons cm~^ s^^, hydrogen density 
TT-H = 10^ cm~^, and hydrogen column density A^h = 10^^ cm~^). The flux of the blackbody is 
varied. In Figure O we plot the Korista continuum and the soft blackbody with several values 
of flux. We understand that the net continuum shown in Figure [6] is not representative of any 
real continuum; this exercise is simply meant to demonstrate the effect of the continuum at 
^ 1640 A, the wavelength of the UV3j transitions (Table [1]), upon the X-ray line ratios. In 
Figure [TJ we plot the i?-ratio as a function of the flux of the second continuum component, 
holding everything else constant. The ratio changes by several orders of magnitude. The 
blackbody only contributes to the net continuum for energies less than about 10-30 eV. This 
is much less than the 165 eV needed to ionize ions in the n = 2 shell of 0^+ and the 560 eV 
needed to ionize ions in the ground state. The blackbody continuum is, however, hard 
enough to directly pump the UV3j lines. This pumping changes the relative populations of 
the 2s ^5* and 2p ^Pj levels, which in turn changes the i?-ratio. 

For large values of the flux of the blackbody, the ratio goes to a well-defined limit as 
the relative populations of 2s and 2p ^Pj levels become populated according to statistical 
weights. In this limit, ^2^35 ^ n2p3p-^, and since A/ ~ Aj, Equation [S] reduces to 

A2,3c_i„lc 1.0 X 10"^ S^^ f- , , 

jl=^lAJl± = - = 1.8x10"^ 18 

^2p3p,_i,i5 5.5xl08s-i ^ ^ 

This limit is also being approached in the extreme upper portion of Figure HJ in calculations 
with very high ionizing flux for similar reasons. 

We also plot the G-ratio in Figure [71 Pumping in the singlet levels by the continuum 
at UVl is much less effective than in the triplet case because the UVl line strength is 
many orders of magnitude less than the resonance line strength, especially considering the 
fact that the resonance line is already being pumped in the low flux limit. Because the 
resonance intensity (the G-ratio denominator) is not strongly affected by the UVl pumping. 
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and because the UV3j pumping decreases one term of the G-ratio numerator while increasing 
the other, the G-ratio is much less affected than the i?-ratio. The G-ratio will finally go to 
a high-flux limit at fluxes roughly three orders of magniude higher than those plotted in 
Figure [3 

4. Conclusions 

The i?-ratio measures the flux in the continuum in a photoionized plasma. It is a good 
density indicator only for a known ionization fraction. The ratio is strongly dependent 
upon continuum pumping of the UV3j lines between n = 2 levels. Because of this 
dependence, X-ray observations can be used to deduce the UV radiation field striking 
the gas where the lines form. This is important for AGN since several components 
of the continuum are beamed, and the continuum we observe may not reflect the 
continuum striking clouds. 

Optical depth in the intercombination line causes the high observed values of the i?-ratio. 
This was previously unexplained. Caution must be taken in using i? as a density 
indicator in conditions where the intercombination line may be optically thick. 

While the G-ratio does track ionization fraction well, the column density must be con- 
strained before using G as a ionization indication, and should not be used as a tem- 
perature indicator at all in a photoionized plasma. 

The resonance line is directly pumped by the incident continuum for small column densities. 
This affects the L and G ratios by large amounts. 

Absorption of UV3j lines may occur in some environments with large column densities. 

For some helium-like species, the optical depth of UV3j lines (or column density of the 
ion) may be determined from X-ray observations of that ion's intercombination lines, 
provided the intercombination doublet can be resolved. The column density could be 
deduced from this. 



We thank Jack Baldwin and Andy Fabian for helpful discussions, and the referee, Manuel 
Bautista, for an insightful review. We also thank NSF for support through AST-0607028 
and NASA for support through NNG05GG04G. 
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A. Data for Helium-like Ions 

Here we discuss the sources and methods of our theoretical models. Where details are 
omitted, the model ion is treated in the same way as the model helium atom (Porter et al. 
2005, Bauman et al. 2005). 

A.l. Energies 

For level energies up to and including n = 5, we use the energies from version 5.0 of 
Chianti (Dere et al. 1997, Landi et al. 2006). Where Chianti does not provide an energy 
we calculate energies using fits to quantum defects and Equation 1 of Bauman et al. (2005), 
modified to include the multiplicative factor {Z — 1)^ in the numerator. This multiplicative 
factor represents the approximate scaling with Z of ionization energy of an ion in the ground 
state. The one is subtracted from the nuclear charge because the charge "seen" by an excited 
electron is partially screened by the charge of the second electron (in the ground state) . In 
this formalism, the quantum defect is nearly independent of Z. 

A. 2. Collisional Data 

For electron impact collisions from ground to n = 2 levels and between the n — 2 
levels, we use simple fits to the effective collision strengths of Zhang & Sampson (1987). The 
most accurate data for these transitions (for oxygen) are from Delahaye & Pradhan (2002) 
who say their results generally agree well with the Zhang & Sampson results. For Stark (l- 
mixing) collisions, we use the method of Seaton (1962) for / < 2 and Vrinceanu & Flannery 
(2001) for greater /. We apply the method of virtual quanta (Jackson 1999) to calculate 
/-mixing collision rates due to proton and singly- ionized helium impact. (These collisions 
are treated in the same way as in the case of helium, described in Porter et al. 2005, where 
they are discussed in some detail.) For n-changing collisions, we use the method of Vriens 
& Smeets (1980), Equations 14-16. We apply the method of virtual quanta to the Vriens 
& Smeets method in order to calculate n-changing collision rates due to proton and singly- 
ionized helium impact. (It is worth noting, however, that electron impact excitation usually 
completely dominates over proton impact for An > transitions.) Collisional de-excitation 
is included via detailed balancing. For collisional ionization, we take the minimum positive 
result from the hydrogenic methods of Allen (1973) and Sampson & Zhang (1988). Proton 
impact excitations are negligible for oxygen but increasingly important for increasing nuclear 
charge (i.e., our rates generally agree with the trend reported by Blaha [1971]). As a basic 



- 15 - 



test of our results in a collisional plasma, we note that our O VII i?-ratio curve at 500,000 K 
has very nearly the same density dependence as is seen in Figure 3 of Bautista & Kallman 
(2000) and Figure 8 of Porquet & Dubau (2000). 

A. 3. Recombination Coefficients 

We add the results of equation 15 of Seaton (1959) for 101 < n < 1000 to the direct 
recombination into the highest collapsed level at n = 100. The recombination coefficients 
into the terms with n < 100 are calculated using the Milne relation. For n < 10 we use fits 
to TOPbase (Fernley et al. 1987; Cunto et al. 1993) photoionization cross-sections for the 
least hydrogenic S terms. We use hydrogenic cross-sections for the rest. Photoionization 
from all levels is included, as is induced recombination via detailed balancing. 

State-specific dielectronic recombination coefficients for n < 8 and / < 4 are interpolated 
from Badnell (2006). Rates to levels with n > 8 are assumed to follow a law, while rates 
to levels with I > 4 are neglected. We make the simple approximation that all doubly-excited 
levels have the same energy, equal to the energy of the n = 2 shell of the H-like stage, so 
that satellite line intensities are proportional to the state-specific dielectronic recombination 
coefficients. We intend to improve this treatment in future work, but the current approxi- 
mation is sufficient for photoionized plasma calculations. We note that Bautista & Kallman 
(2000) warn that satellite lines lie close in wavelength to the forbidden and intercombination 
lines and will enhance their apparent intensities in low and medium resolution spectra. 

A. 4. Transition Probabilities 

A.4-1- Allowed Transitions 

We use the transition probabilities of Johnson et al. (2002) where available. For n = 2 
to n = 2 transitions not available in Johnson et al, we use simple fits to data from the NIST 
Atomic Spectra Database, version 3.oM for 2p%-2s^S and TOPbase for 2p^P-2s^S. We 
also use TOPbase for same-n transitions with n> 2. Then, for transitions with either initial 
or final / > 2, we use the hydrogenic formula given by Drake (1996). Next, transitions to 
ground, 2s ^S, 2s ^S, 3s ^S, or 4s are calculated via extrapolation of Johnson et al. values 
to higher initial n. All other allowed transitions are calculated through the use of Drake's 



•^see ,http://physics. nist.gov/PhysRefData/ASD/index.htinl 
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semi-classical quantum-defect routine. 



A.4-2. Forbidden Transitions 

The transition probabilities for the transitions 2s '^5' — Is are calculated from a fit to 
Lin et al. (1977) values. Two-photon transition probabilities are taken from Derevianko & 
Johnson (1997). The distribution of photon energies are from Johnson (2002). The sources 
of inter combination transition probabilities are as follows: upper level 2p^Pi: fits to Johnson 
et al. (2002) for Z < 18, and Lin et al. (1977) otherwise; and upper level 2p^P2: fits to Lin 
et al. (1977). Transition probabilites for transitions np^P — n's^S and np^P — n's with 
n ^ n' are taken from fits to (or extrapolation of) Johnson et al. values. Finally, transition 
probabilities for the transition 2p ^P — 2s are taken from fits to values in Savukov et al. 
(2003). 



B. Absorption of Intercombination and UV3j Lines 

It is not generally known that absorption of the intercombination and UVSj lines of 
some helium-like ions may be possible in some environments. Here we present arguments in 
this regard and derive column densities necessary to detect the absorption of these lines. 

The line-absorption coefficient of a line with lower and upper levels / and u and central 
frequency uq, neglecting radiative damping, is 



A — 



2 



cm% (Bl) 



where e is the electronic charge, me is the electronic mass, c is the speed of light, fiu is the 
absorption oscillator strength of the line, and Aiy^ is the Doppler width [Hz] (which consists 
of thermal and turbulent components added in quadrature). Equation IBll reduces, at line 
center, to 

«:o = 0.0150%^ [cm^]. (B2) 
The optical depth at line center is 

To = y kq ni{Z,T) dl [dimensionless] , (B3) 

where ni{Z, T) is the population of level I and the integral is over the length of the absorbing 
region. 
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For the 2p ^Pj — 2s UV3j transitions of helium-like ions (with wavelengths given in 
Tabled]), we substitute Equation [5] into equation IB3I and, assuming an isothermal absorbing 
region, obtain 

To^Ko f'^'^^'^X, [n^n^.dl. (B4) 



^2s3s'_isis(Z) 

If we further assume that the electron density does not change too much over the length of 
the absorbing region, then 

Tq^Kq-. 77^ l\z+, lB5j 

where Nz+ is the column density (cm~^) of the hydrogen-like ion. 

For oxygen at 500, 000 K, we find al{,ig{8, 10^-'^ K) ^ l.lx 10"^^ cm^ s"^ and, neglecting 
turbulence, Au^ = 1-4 x lO^^Hz, so that tq ~ 5.7 x 10^^° Ue No'r+ for the strongest line of 
the triplet. For appreciable absorption, letting tq = 0.1, we need Ue No7+ > 2 x 10^^ cm~^. 
At Tie = 10^° cm~'^ one would need No'r+ ~ 2 x 10^*^ cm^^ to detect absorption. That column 
density is one to two orders of magnitude larger than is typically reported but may exist 
in some environments. A larger volume density would allow for a smaller column density. 
Absorption is thus likely to be detectable in some environments. The likelihood of detection 
decreases, however, with increasing nuclear charge for two reasons. First, the transition 
probability in the denominator of Equation IB5I increases much faster than any other factors 
change. Second, astrophysical abundances tend to decrease dramatically for elements heavier 
than oxygen. Absorption of the UV3j lines will be considered in a future paper. 

For the 2p^Pj — Is^S transitions of helium-like ions, we can make the approximation 
that all ions are in the ground state and obtain 

To ^ Ko Nz, (B6) 

where Nz is the column density (cm~^) of the helium- like ion. Again considering oxygen at 
500, 000 K, we find Au^ = 1.1 x lO^^Hz, so that to^2.1x IQ-^^ No6+ for the strongest line 
of the doublet, where we have used fis^s,2p»Pi = 1-2 x 10^"^. Again letting tq = 0.1, we need 
-^06+ > 5 X 10^^ cm~^. If only one-tenth of the oxygen in the column is helium-like oxygen, 
and the oxygen abundance was comparable to solar, the total hydrogen column density 
necessary to see significant absorption of the O VII intercombination line would be about 
10^^ cm~^. Absorption is thus likely to be detectable in some environments, which has been 
stated by Godet et al. (2004), and may have been detected (see Kinkhabwala et al. 2002). 
Including turbulence will reduce the optical depth and increase the necessary optical depth 
to observe absorption (although the effects will be strongly dependent upon line profile). 
The likelihood of detection decreases with increasing nuclear charge, however, because the 
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intercombination rate increasingly dominates the UV3j rates (so that the intercombination 
lines effectively become the only decay mode from 2p ^Pj) , and because abundances tend to 
decrease dramatically for elements heavier than oxygen. 
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Table 1. Wavelengths of the UV, Intercombination, Forbidden, and Resonance Transitions of He-hke Ions. 



Wavelength (A) 
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UV3i 


UV32 


UVl 




f 




kJ V 1111^^1 


z 
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2270.8 
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N 


7 
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1907.3 


1896.7 


2896.4 


29.084 


29.534 


28.787 


O 


8 


1639.9 


1638.3 


1623.6 


2449.0 


21.807 


22.101 


21.602 
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1417.2 


1414.2 


1395.5 


2136.1 


16.947 


17.153 


16.807 


Ne 
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1277.7 


1272.8 


1248.3 


1856.0 


13.553 


13.699 


13.447 


Na 
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1149.2 


1142.3 


1111.8 


1646.9 


11.083 


11.192 


11.003 


Mg 


12 


1043.3 
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997.46 


1474.2 


9.2312 


9.3143 


9.1688 


Al 


13 


954.33 


943.16 


899.67 


1327.7 
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Si 


14 


878.65 


865.14 


814.69 
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Ca 
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Fe 


26 
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271.16 
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1.8682 
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Fig. 1. — Grotrian diagram of the n = 1 and n = 2 levels of helium-like oxygen. Some 
transitions are not shown. Relative energies are not drawn to scale. The energy order is 
different for different ions of the sequence. 
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log (A/|-|/cnn ) 

Fig. 2. — R = li/Ii, L = Ir/Ii, and G = (If + /i)//r as a function of hydrogen column density. 
See text for details of model. 




Log /iH {cm"^) 

Fig. 3. — Ionization fraction of oxygen in the (left panel) O^"*" and (right panel) O^^ stages. 
Contours are labeled as the log of the fraction. The column density is A^h = 10^^ cm~^ 
in both panels. The case with A^h = 10^^ cm~^ is very similar except that the lower 10~^^ 
contour is shifted up nearly a decade for 0^+ and about half a decade for 0^+. Wavy features 
in this and following contour plots are artifacts in the plotting program due to a sparse grid. 
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Log (cm"?) 

Fig. 4. — Ratio R — If/Ii as a, function of hydrogen density and ionizing flux. Left panel 
is for a small column density (A^h = 10^^ cm~^); right panel is for a large column density 
{Nil = 10^^ cm~^). The bold diagonal lines in each panel show constant ionization parameter. 
See text for details of model. 
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Fig. 5. — Same as Figure H] except ratio is G = (If + /i)//r- Note that the contours have a 
hnear scale in the left panel and a logarithmic scale in the right panel. When the ionization 
fractions (Figure [3]) peak, the G-ratio is sensitive to ionization parameter and column density. 
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Fig. 6. — A typical AGN continuum is shown as the sohd hne. The soft blackbody, which 
we vary to show its effect on the X-ray hues, is also shown for several different values of flux. 
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Fig. 7. — i?-ratio and G-ratio as functions of the flux of the soft component shown in Figure [61 
The low flux limit is the canonical value. The i?-ratio high flux limit is where the relative 
populations of the 2s and 2p^Pi levels are statistically weighted. The corresponding high 
flux limit in the G-ratio has not been reached because the UVl line strength is several orders 
of magnitude less than the resonance line strength. 



